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Abstract 

We define a current-conserving approximation for the local conductivity ten- 
sor of a disordered system which includes the effects of weak localization. 
Using this approximation we show that the weak localization effect in conduc- 
tance is not obtained simply from the diagram corresponding to the coherent 
back-scattering peak observed in optical experiments. Other diagrams con- 
tribute to the effect at the same order and decrease its value. These diagrams 
appear to have no semiclassical analogues, a fact which may have implications 
for the semiclassical theory of chaotic systems. The effects of discrete sym- 
metries on weak localization in disordered conductors is evaluated and and 
compared to results from chaotic scatterers. 
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I. INTRODUCTION 



The weak localization (WL) effect has played a seminal role in both theoretical and ex- 
perimental work on quantum transport in disordered conductors. On the theoretical side it 
enters the scaling theory of localization in a fundamental manner and on the exper- 
imental side it provides the best measurement of the phase- coherence length for electrons 
at low temperature ||. Finally, through the normal-metal Aharonov-Bohm effect, it pro- 
vides a very direct measurement of quantum interference of normal electrons on the scale 
of this phase-coherence length ||. The quantitative theory of weak-localization is based 
on impurity-averaged perturbation theory in the small parameter (kfl) -1 (kf is the fermi 
wavevector and I is the elastic mean free path) [|J. Many experiments over the past decade 
have tested this theory in detail and found excellent agreement with the measured low- 
temperature variation of the conductance with temperature and magnetic field. 

The weak localization effect arises due to time-reversal (TR) symmetry. It is often de- 
scribed by noting that a diffusing electron wave has an enhanced probability of returning 
to its starting point (compared to a classical particle) due to the constructive interference 
of time-reversed pairs of trajectories forming closed loops. This enhancement of the return 
probability then leads to reduced conductance with respect to the classical Drude value. 
Breaking time-reversal symmetry eliminates this interference causing a positive magneto- 
conductance (in the simplest case in which spin-orbit scattering is negligible ||). Although 
WL is an interference effect and hence clearly non-classical, it is often described as a semi- 
classical effect 0,[5]|| which can be understood simply by adding the amplitudes for motion 
along classical trajectories with appropriate phases, and then squaring the amplitude. Since 
recent theoretical advances in semi-classical quantum theory have shown quantitative suc- 
cess in other fields, e.g. atomic physics 0; it is of some interest to see if weak localization 
can be quantitatively obtained from a summation over classical paths. Recently two of the 
authors addressed this problem in a treatment of WL in ballistic conductors and numeri- 
cally found contributions to the WL effect which they could not obtain from a semiclassical 
calculation f|,[J. It was noted || that a similar contribution apparently entered in the WL 
theory of disordered conductors, but had not been explicitly discussed or calculated in the 
literature. We present such a calculation below and discuss its implications. Reference || 
also discovered a sensitivity of WL in ballistic conductors to discrete spatial symmetries; we 
evaluate this sensitivity for disordered conductors below. 

Well before the discovery of weak localization in disordered conductors the increase in 
backscattering probability due to time-reversal symmetry had been noted both in optics 



P0| . |Tl|| and in nuclear scattering ||12|| . In the former this effect is referred to as coherent 
backscattering (CB) and in the latter case as the elastic enhancement factor. In optics one 
clearly observes a peak approaching a factor of two in height in the backwards direction 
[r5HH)| (if the polarization is preserved). Recently the CB effect in optics has been hypoth- 
esized [[UJ to explain the surge in brightness of astronomical bodies when in opposition, and 
this hypothesis is believed to have significant implications for planetary science. Although 
known for some time, the detailed study of CB in optics was stimulated by WL theory and 
the two terms (weak localization and coherent backscattering) are often used interchange- 
ably. In this work we emphasize that the two effects are not equivalent: diagrams which 
make a negligible contribution to the differential scattering cross-section make a contribution 
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to the weak localization conductance of the same order and of opposite sign as the diagram 
which gives rise to the CB peak. 



II. CURRENT CONSERVATION FOR EXACT CONDUCTIVITY AND 

CONDUCTANCE 

In linear response theory the quantum expectation value of the local DC current density 
induced by an electric field E(r') is 

< J( r ) >= J drV(r, r') • E(r') (1) 

where a is the Kubo local conductivity tensor, which may be expressed in terms of the exact 
eigenfunctions and eigenenergies for a given disorder configuration [|T7|-|T9||. The general 
expression for g_ in the presence of a magnetic field, B, may be divided into a term which 
only depends on states at the fermi surface (at T=0) and is symmetric in B, and a term which 



depends on all states below the fermi surface and is anti-symmetric in B ||19|| . The weak 
localization effect relates to average properties of conductivity and conductance which are 
symmetric in B, hence the anti-symmetric term may be neglected. With this simplification 
we may write 

- (r ' r/) = " ii^ AG(£/ ' r ' r ' } ^ * v AG ( £ f> r '' r )' ( 2 ) 

where AG(E,r,r') = G + (r,r') — G~(r, r'), G ± is the advanced (retarded) Green functions 
for the single particle Schrodinger equation, and D represents the antisymmetric gauge- 
invariant derivative: f(x)Dg(x) = f(x)[Dg(x)] — g(x)[Df(x)], D = V — (ei/hc)A(r). Since 
WL occurs at relatively weak magnetic fields it turns out to be sufficient to make the further 
approximation D « V which we do henceforth. With this set of approximations DC current 



conservation for the local current density, V- < J(r) >= 0, is equivalent [19] to the condition 
V-a = V'-a = 0. 

Transport experiments do not directly measure the local conductivity tensor; instead 
they measure conductance or resistance. In the simplest case of a two-probe measurement, 
in which the voltage drop is measured between the current source and sink, the conductance, 
g, is just the inverse of the resistance and is related to the local conductivity by integration 
over the cross-sections at the interface between the sample and the leads [|19|j : 

g = dSi / dS 2 hi ■ a(r, r') • n 2 (3) 



where ni,n 2 are unit normal vectors to the cross-sections of the two leads. 

There are two different forms of Eq. (3) often used in the literature; the first reexpresses 
conductance in terms of scattering coefficients, and the second reexpresses conductance as an 
average over the sample volume. In the first case, the cross-sections are taken just outside 
the sample in order to express g_ in terms of the exact asymptotic wavefunctions which 
involve the transmission and reflection amplitudes. Straightforward manipulations worked 
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out by Economou and Soukolis |L7[] and by Fisher and Lee Jl8| then show that Eq. (3) is 
equivalent to the two-probe Landauer formula: 



S = tE3« (4) 

where is the transmission coefficient between the propagating states i,j in the leads 
(assumed identical) and there are N c such states at Ef. Below we will discuss extensively 
the contribution to the WL effect of the transmission coefficients which are diagonal or 
off-diagonal in the mode indices 

In the second case one uses the relation V • a — V • g_ — to show that g is independent 
of the location of the cross-section in Eq. (3) and then integrates over both x, x' to obtain 
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J drdr'AG(e f , r, r')V x V x ,AG(e f , r', r), (5) 



where the integrations are over the entire sample. This volume-averaged form for the con- 
ductance has been used in most previous studies of WL; as we will see it leads to a much 
simpler set of diagrams for the average conductance than does Eq. (4). We note that the 
volume-averaged form is only valid for two-probe conductance, and treatments of four-probe 
resistance measurements must rely on Eqs. (3)and (4) [|I9|-|2T||. This is particularly impor- 



tant for mesoscopic samples, where the lead geometry may have significant influence on the 
measurements. 



III. CURRENT-CONSERVING APPROXIMATIONS FOR < a > 

A. Impurity- Averaged Perturbation Theory 

Disordered conductors are typically characterized by their statistical properties when 
averaged over an ensemble of impurity configurations. As noted above, the WL effect appears 
in the average conductance; the primary analytic tool for calculating this effect is impurity- 
averaged perturbation theory in the small parameter (kfl)^ 1 . This perturbation theory is 
equivalent to that for electrons with a static interaction, except for the absence of closed 
fermion loops (which is explained naturally in the replica formulation of the problem). As in 
any diagrammatic expansion, it is crucial to evaluate a set of diagrams at each order which 
maintains the conservation laws present in the exact theory. In this case the important 
conserved quantity is current, and we seek a set of diagrams for < a(r, r') > which both 
describes the WL effect and satisfies 

V- < a>= V- < a >= 0, (6) 

where henceforth both the angle brackets and overbars will denote the average over impurity 
configurations. 

We shall adopt the standard white-noise model for the statistics of the impurity potential: 
< V(r) >= 0, < V(r)V(r') >= Cj-u 2 5(r — r'), where q is the density of impurities and all 
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higher cumulants are assumed zero. The usual approximation (see e.g. Ref. p8| ) for the 
average Green function is the self-consistent Born approximation (SCBA) which is shown 
diagrammatically in Fig. 1; it may be expressed for our case by the self-consistent equation: 

[E - Ci u 2 G(r, r) + ^]G(r, r') = 5(r - r'). (7) 

Although we will consider finite disordered samples connected to ordered leads (and this will 
be important for certain aspects of the calculation) the violation of translational symmetry 
at the boundaries is negligible for the average Green function. This is first because G decays 
exponentially for |r — r'| > I and is insensitive to the boundaries over most of the sample, 
and secondly because even at the boundaries corrections to the average Green function due 
to the interface only appear at next order in (kfl)^ 1 . Since the boundary effects are found 
to be negligible in Eq. (7) we may treat the system as translationally invariant on average, 
hence G(r, r, E) is just a complex function of energy, independent of r. For weak magnetic 
field even the energy-dependence is negligible and one finds CiU 2 G (r, r) = =)=i/2r, where r 
is the elastic scattering time. It then follows from Eq. (7) that 

G ± {v,v',E) = G±{r,T , ,E±i/<lT), (8) 

where G is the Green function of the problem in the absence of disorder. Substituting this 
relation into the known form of Gq leads to the exponential decay with I mentioned above. 



B. Current Conservation in the Ladder Approximation 

At leading order in (kfl)" 1 it is well-known [2~2~| that the SCBA for the average Green 
function will yield a current-conserving approximation for the conductivity if the velocity 
vertices are dressed by the infinite sum of ladder diagrams (see Fig. 2). This ladder sum 
L(r, r') satisfies the integral equation: 

L(r, r') = 5(r - r') + J d ri L (r, n)£(n, r') (9) 

where the kernel L (r, r') = CiU 2 \G + (r, r') | 2 . 

The proofs of current conservation in various approximations to < a(r, r') > rely on 
properties of the velocity vertex under the divergence operation. The most general (bare) 
velocity vertex which can appear has the form 

V (r, r 1; r 2 ) = G + (r 1; r)^G-(r, r 2 ) (10) 

where r 1; r 2 are normally integrated over in evaluating the diagram and V is the anti- 
symmetric gradient operator. To prove current conservation for any set of diagrams we 
must take the divergence of each diagram in the chosen set and show that the resulting sum 
vanishes. Since the r, r' dependence in each diagram resides solely in the vertices we need 
only evaluate the divergence of the relevant vertex and then convolve this with the rest of 
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the diagram. Following Hershfield p2| , we can regard the operation of taking the divergence 
as generating a new set of diagrams with the velocity vertex replaced by a "divergence 
vertex" which simply designates the point at which the divergence was taken (we denote 
this divergence vertex by a black dot in the figures) ; the resulting set of diagrams must sum 
to zero if the approximation conserves current. 

The divergence of the vertex Vq may be evaluated using Eq. (7) to substitute for the 
factors V 2 G/2m which arises after differentiation. 

V ■ V {t, n, r 2 ) = -G + ( ri , r)G"(r, r 2 ) - iG~(r, r 2 )5(r - n) + zG+(r a , r)5(r - r 2 ). (11) 
r 

The diagrams generated by this procedure are shown in Fig. 3. We see that in general 
taking the divergence of a bare vertex in a conductivity diagram generates three diagrams: 
the first and second diagrams have one of the Green functions deleted and the divergence 
vertex moved to r 1; r 2 respectively, and the third is identical to the original conductivity 
diagram except that the divergence vertex replaces the velocity vertex and is multiplied by 
a factor 1/r, 

Now consider dressing the bare velocity vertex with an impurity line as shown in Fig. 4a. 
In this case the bare vertex connects to the same point, r ; if we set ri = r 2 = r in Eq. (11) 
the second and third terms combine using the identity — iCiU 2 [G (r, r) — G (r, r)] = —1/r 
and yield a single diagram with the impurity line removed as shown in Fig. 4a. A similar 
result occurs when the divergence is taken of a bare vertex dressed by two impurity lines. It 
is evident from the figure that the first diagram generated from the singly dressed vertex is 
the negative of the second diagram generated from the doubly dressed vertex and will cancel 
if these two diagrams are summed. It is very useful at this point to define the dressed vertex 
V as the vertex V dressed with a full ladder sum (note that the bare vertex appears as the 
first term of V). Because of the cancellation of successive diagrams generated by taking the 
divergence of a velocity vertex connected to a ladder we obtain the simple result 

V • V(r, n, r 2 ) = -iG"(r, r 2 )6(r - n) + zG + (n, r)5(r - r 2 ). (12) 

This result is shown diagrammatically in Fig. 4b. The divergence operation on the dressed 
vertex generates only two diagrams, one with the point r shifted to the point ri, the other 
with it shifted to the point r 2 , with the analyticity of the Green function reversed, and with 
opposite sign. These two diagrams will have the same number of impurity lines. 

With this identity it is easy to confirm current conservation in the ladder approxima- 
tion for < g_ > (this approximation is represented diagrammatically in Fig. 5). Taking the 
divergence at the dressed vertex generates two diagrams with a single Green function multi- 
plied by 8(r — r') attached to a bare velocity vertex. The antisymmetric derivative of these 
diagrams vanishes due to the symmetry of the Green functions under interchange of r, r'. 
Thus, as claimed, the divergence of < g_ > in this approximation is zero. It is worth noting 
that this argument establishing current conservation remains valid even at the interface of 
between the sample and the leads. 
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C. Current Conservation Including Weak Localization 



The ladder approximation to < a > discussed in the last section does not, however, 
describe the WL effect: the WL effect only appears at next order in (kfl)~ l . We have 
not found in the literature an explicitly current- conserving set of diagrams for < a > which 
upon integration yields the WL contribution to the conductance. However, it is known || 
that the volume-averaged expression for g yields the full effect when one evaluates only the 
diagram with one crossed ladder (Fig. 2b) without vertex corrections; this implies that all 
other volume-averaged diagrams at this order either vanish or sum to zero. Thus the WL 
effect in the total conductance may be calculated without determining an explicit current- 
conserving set for < g_ >. However in order to evaluate individual transmission or reflection 
coefficients < >,< > and determine their separate contributions to < g >, the 
identification of such a set is necessary. A current-conserving approximation which includes 
the WL effect is represented by the diagrams shown in Fig. 6. This approximation will be 
shown to conserve current exactly (if the SCBA is used for the average Green function) by 
a generalization of the type of argument p^] used for the simple ladder approximation to 
< g_ > and for conductance fluctuations in multi-lead structures. The set given includes 
diagrams of higher order in (kfl)" 1 than needed for the WL effect; however, inclusion of 
the higher-order diagrams actually expedites the proof. (We will not evaluate all diagrams 
explicitly anyway). In these diagrams the crossed insertion represents the time-reversed or 
crossed ladder shown in Fig. 2b, while the number on each insertion represents the number 
of impurity scatterings. We will show that the divergence of this set for a fixed total number 
of impurity lines is zero; summing over these sets for all numbers of impurity lines in the 
ladder yields a complete set of current- conserving diagrams containing the crossed ladder 
sum (and hence the WL effect). 

The proof that the set of diagrams shown in Fig. 6a,b satisfy V- < a >= for each n is 
as follows. The set shown consists of a crossed ladder with n scatterings and all diagrams 
containing n — j scatterings in the crossed ladder and j crossed lines dressing the crossed 
ladder on one of the two Green's functions, for j = 1, 2 ... n. (For j > 1 we have included 
diagrams of lower order in (kfl)^ 1 .) It is helpful to redraw these diagrams as shown in 
the figures; in this way it is clear that this set is generated simply by moving one of the 
vertices V around the figure, passing through the impurity lines one by one. Without loss 
of generality, we may take the divergence at the "moving" vertex. As just shown, taking 
the divergence of each diagram generates two diagrams, one in which the vertex is moved 
to the impurity line immediately to the left, the other in which the sign is reversed and the 
vertex is moved to the impurity line immediately to the right (and in both cases the ladder 
dressing is removed). This alternation in sign leads to a cancellation between the "right" 
diagram generated from the conductivity diagram with j = m and the the "left" diagram 
generated by the one with j = m + 1. This cancellation precedes around the figure until 
the only remaining divergence diagram is the one where the divergence vertex has moved 
all the way around the crossed ladder and lies adjacent to the remaining dressed velocity 
vertex attached at r'. This diagram may be exactly cancelled by adding to the conductivity 
the set of diagrams shown in Fig. 6b consisting of a crossed ladder self-energy insertion 
to one of the Green functions dressed by an ordinary ladder. The divergence of this new 
set cancels internally except for the diagram shown in Fig. 6b. This diagram is exactly 
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the negative of the remaining divergence diagram from the first set, except that it doesn't 
contain the bare vertex. Thus the sum of the divergence of both sets cancels exactly, except 
for the diagram with the bare velocity vertex evaluated at r = r'. However we noted above 
that such a diagram vanishes when the antisymmetric derivative with respect to r' is taken. 
Thus we have shown that the sum of the two sets does indeed yield a current-conserving 
approximation to < a >. 

We note that this operation of "moving a vertex around the figure" is a quite general 
way of constructing a current-conserving set of diagrams from a given diagram. In fact, 
the additional ladder dressings on the crossed ladder that were needed to complete the 
cancellation can be regarded as the result of moving the vertex through the ladder dressing 
on the other vertex V (the one attached to the point r' in Fig. 6a. 



D. Limitations of Diffuson and Cooperon Approximations 

The exact ladder sum would be obtained by solving the integral equation (9), however 
in most applications to disordered conductors one is interested in the solution in the limit 
of long wavelengths, so one assumes that L(r, r') is slowly- varying compared to L and one 
expands L(r 1 , r') for r x « r to obtain a differential equation for L(r, r') which has the form 
of a diffusion equation [2S|. For the simplest case of zero frequency and temperature this 
takes the form 

-DrV 2 d(r,r') = 5(r-r'), (13) 

where D is the elastic diffusion constant, l 2 /dr, and we denote the long wavelength ap- 
proximation to the ladder sum by d(r, r'). We will henceforth refer to the exact solution of 
Eq. (9) as the ladder approximation or ladder sum while the solution to Eq. (13) will be 
referred to as the diffuson approximation (the terms are often used interchangeably in the 
literature). In the absence of a magnetic field the crossed ladder sum of Fig. 2b satisfies 
the same integral equation as the ladder sum, and thus the long-wavelength approximation, 
which we denote as the cooperon, satisfies the same diffusion equation. However, in contrast 
to the diffuson, the cooperon is sensitive to weak magnetic fields, which enter the diffusion 
equation by the minimal substitution for particles of charge 2e, iV — > zV — 2eA/hc P,|28" 



The appearance of the magnetic field in this way leads to the suppression of weak localization 
by a magnetic field through the elimination of the diffusion pole. 

Although the diffuson/cooperon approximation to the ladder /crossed-ladder sums is ad- 
equate for evaluating the conductivity tensor for widely separated points r, r', it causes a 
breakdown of local current conservation. With appropriate care this breakdown can be 
shown not to affect the value of the average transmission coefficients and hence of the con- 
ductance, essentially because one is only evaluating the diffuson at widely separated points. 
However the diffuson/cooperon approximation is not adequate to calculate the average re- 
flection coefficients (a fact which has been known in the optics literature for some time |[23|| ). 
The reason for this is that a major contribution to reflection involves short trajectories of 
only a few mean free paths in length, for which one needs an accurate expression for L(r, r') 
on the scale of I. But the exact short wavelength behavior of the ladder L(r, r') is known (for 
an infinite medium) [El] and differs substantially from that of the diffuson approximation. 
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As a result, although the diffuson approximation allows a correct evaluation of < T > it 
fails for < R > (the average total reflection) and hence < T > + < R >^ N c , and the 
local violation of current conservation shows up as a global violation of unitarity. We now 
describe this problem in some detail. 

First we discuss the average total transmission, < T >, which is just < g > /(e 2 /h) 
by Eq. (4). If one starts from the volume-averaged expression for g it is well-known that 
the Drude-Sommerfeld average conductance is correctly obtained from the leading order 
diagram in (kfl)^ 1 within the white-noise model, i.e. one finds < g >= a A/L where 
o" = ne 2 r/m and A, L are the cross-sectional area and length of the conductor. It is easy 
to show that this conductance corresponds to an average transmission < T >= (tt/2)NJ/L 
(in two dimensions). In the volume-averaged approach, however, all the ladder diagrams for 
the conductance vanish from isotropy so the result is insensitive to the failure of the diffuson 
approximation at short distances. 

In contrast, when evaluating conductance from the transmission coefficients the ladder 
diagrams do not vanish since the conductivity tensor is only integrated over a surface. We 
have not found in the literature confirmation that calculating < T >= YliJ < Ty > in 
the diffuson approximation yields the correct result as well. We now confirm that it does, 
although only after improving the standard approximation for the diffuson boundary condi- 
tions. In volume-averaged calculations the diffuson is taken to vanish on the surface of the 
conductor connected to the leads (this is supposed to correspond to the boundary condition 
of a very large metallic contact). However if this approximation is made and if the point r 
is chosen at the surface, then < <x(r, r') > will not conserve current even for r' far in the 
bulk. It can be shown however that if the diffuson is set to zero a distance of order the mean 



free path outside the interface |25j this property is restored for the total flux crossing the 
interface. In two dimensions the correct distance is nl/4 outside the sample. This modified 
boundary condition is required to get the correct average transmission coefficients; to get 
correct reflection coefficients even this approximation is inadequate and a new approximation 
is needed which is discussed below. 

The average transmission and reflection coefficients to leading order are given by the 
diagram consisting of a single diffuson between two current vertices. To examine individual 
Tij,Rij we must evaluate these diagrams for a given modes i,j which requires integrating 
over the transverse directions (y, y') with the mode wavefunctions while staying in real-space 
for the current direction (x,x'). This is the method we use henceforth in the paper when 
evaluating diagrams explicitly; all calculations are done in two dimensions for convenience. 
Upon transforming to mode space, the differential operator at the current vertices V/2mi 
simply is replaced by the longitudinal velocity, 

< Tij >= ViVj < G + (x = 0, i; x = L x , j)G~(x = L x , j; x = 0, i) > (14) 

where v j, Vj are velocities for the relevant modes. We need the 1PGF which becomes 

_l_ r dp gipx{x—x') 

G (X ' ' )= M ( 2? e/ -p72m + i/2r ^ 
= -i5 id e ip ^ x - x ' Hx - xV{2v i T) / Vj (15b) 
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where p 2 is the sum of squares of transverse and longitudinal momenta and pj = 
~ ( fji < lj = U^/Ly) 2 - Finally, the propagation from a given mode at x = to a 
point r' is given using this Green's function and projecting from a mode to a y-position. We 
obtain 

G + (x, j- r') = -i-jvi\*-x'\-\x-*'\I^T) sin ( g ^')( 2 /L J/ ) 1 / 2 (16) 

V 3 

where the factor of (2/Ly) 1 ^ 2 is added for normalization. In evaluating the diagram for 
< Tij > we note that the phase factors involving x, x' cancel since the advanced and retarded 
Green's functions have opposite phases. The phases coming from y, y' vary on the length 
scale kj 1 , while the diffuson varies on at most the length scale /. So, we can replace the 
y-phase factor arising from the square of sin^y') by its average. One finds: 

< >= Cl u 2 -^ £ X dxdx' Lv~^e- x /^e- x '/^ J dy J dy'd(x, y; x', y')) (17) 

where x is measured from the left-hand side of the sample and x' from the right-hand side. 
The integral of the 2d diffuson over y, y' satisfies a Id diffusion equation which has the simple 
solution: d(x,x') — (x + \l){x' + jl)(2L y /l 2 L x ), where the final factor of L y arises from the 
transverse integration and we have dropped terms lower order in (1/ L x ). The integral above 
can now be integrated by parts; we obtain after some algebra and dropping exponentially 
small terms 

< T 3 >= ^U^- + tt/4)(^ + tt/4). (18) 

7TiV L x Vf Vf 



Summing over the modes i,j using the relation mvj = ^p 2 — (jir / ' L y ) 2 yields the correct 
result 

N ° 7T / 

E < ^j >= N c -j-. (19) 

i,j ^ 

If global current conservation were satisfied in the diffuson approximation we should find 
Y^fj < Rij >= N c (l — nl/2L x ) ks N c . However it is shown in the Appendix that exactly 
the same calculation for < > yields the result Y^ij < Rij >~ 0.73A^ C . As noted above, 
this failure of global current-conservation arises from the fact that reflection coefficients 
require consideration of the conductivity tensor on scales of order the mean free path, and 
the diffuson approximation to the ladder sum is not accurate on this scale. 

E. Compensating for the Diffuson Approximation 

The most elegant resolution to this problem would be to work with the full ladder sum 
or some improved approximation to it, instead of the diffuson. A partial improvement to the 
diffuson approximation is given in the Appendix, however it is not sufficient to restore current 
conservation to this order in (kfl)~ l . Although an exact evaluation of the crossed-ladder 
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sum (which differs from the ladder sum only in the absence of the single scattering diagram), 
does exist [Q in momentum space, it is not clear how to obtain a real space solution that 
can be used in a sample with boundaries. The corresponding real space differential equation 
would be of infinite degree, complicating the specification of boundary values. The analogous 
problem in classical diffusion, essentially the Milne problem, has been solved f2l| , but the 
solution appears too complicated to be of any practical use in this context. 

An alternative approach to restoring current conservation is to find an approximation for 
the current vertex when connected to the ladder sum which restores current conservation 
when the ladder is approximated by the diffuson. This was essentially the approach taken in 
the work of Kane et al. [pO| , pT| , where they treat the bare conductivity bubble as proportional 
to <5(r — r'). However calculation of the total reflection (for which the points r, r' may 
coincide) requires a slightly more careful treatment to avoid evaluating this delta function 
at zero. First we note that current conservation in the ladder approximation to < g_ > 
followed from the cancellation between the term (1/t)\G (r, r')| 2 arising from V ■ Vo(r, r') 
and that arising from the divergence of the bare vertex dressed by the ladder sum. We can 
find a current-conserving approximation when the ladder sum is replaced by the diffuson if 
we can modify Vo in the dressed vertex so as to maintain the result that the divergence of 
the dressed vertex is — (l/r)|G + (r, r')| 2 . The required modification can be determined by 
the following argument. Note that the bare vertex depicted in Fig. 3 with ri = r 2 = r' is a 
vector function of r — r' with zero curl. Hence it may be written as the gradient of a scalar 
function F(r — r'). Therefore the divergence of the dressed vertex becomes 

J dr'dr"V 2 F(r - r')d(r', r")G + (r", ri)G~(r", r 2 ). (20) 

Because F is a function of r — r', integrating by parts twice shifts the Laplacian onto the 
diffuson and allows use of the basic equation V 2 d(r" — r') = — (1/Dt)5(v" — r') to give the 
result 

-j^j dr"F(r-r")G + (r", ri )Glr",r 2 ). (21) 

It follows that current conservation requires F(r — r') = D5(r — r') or equivalently that 
V-F = Vo(r, r') = DV5(r-r'). Of course Vo(r — r') is known explicitly; it is an antisymmetric 
function of range I, not of zero range. Hence as claimed above the diffuson approximation 
does not exactly conserve current. However it is possible to let the mean-free path tend 
to zero in the vertex Green functions such that current conservation is restored; i.e. such 
that the vertex behaves as the gradient of a delta function. We note that for reflection it 
is not possible to let all bare vertices approach the gradient of a delta function because of 
the necessity of evaluating the vertex at r = r'; only the vertex connected to the diffuson is 
taken to be of zero range. 

The desired limit is achieved formally by changing the mean free path in these Green 
functions from / to I' — > 0, rescaling kf to (kf)lfl', and rescaling the Green functions by the 
overall factor l/V. The antisymmetric derivative operator makes the vertex in direction i 
an odd function of (rj — r ■) and an even function of other coordinates, where is the i th 
position coordinate of the derivative operator and r- is the ith position coordinate of the 
other point to which the Green's functions in the vertex are connected. It is straightforward 
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to check that the new vertex vanishes when the function against which it is convoluted 
vanishes to zeroth and first order in any coordinate, and that the vertex in direction i does 
not vanish when convoluted against a linear function of (r* — rQ. This means that the vertex 
acts on a polynomial in the coordinates to select out only terms linear in one coordinate and 
constant in the other two. Therefore the vertex does act as the gradient of a delta function, 
and it is again easy to check that the weight of the delta function is exactly the diffusion 
constant, D, as required for current conservation. A quick way of confirming that the weight 
is correct is to note that for points on opposite sides of the sample the difference between 
using a vertex of range I and the new zero range vertex introduces only exponentially small 
corrections in l/L x . Since we confirmed above that the diffuson approximation with the old 
vertex gives exactly the correct value for < T > to leading order, it follows that the leading 
behavior of the new vertex must be quantitatively correct. 

Using this limiting procedure on the bare vertex allows us to calculate the correct average 
total reflection probability using the diffuson approximation (and employing a real-space 
formulation analogous to Eq. (3)), i.e. we find < R >= N c to leading order in l/L x . 
However this procedure has significant limitations. It only works for evaluation of real-space 
diagrams or diagrams with one vertex in real-space and one in mode space (the one not 
taken to be zero range). We have not found a sensible current-conserving approximation 
incorporating the diffuson approximation which will allow us to evaluate the individual 
< Rij >, so explicit calculations of reflection coefficient diagrams typically contain errors of 
order unity in prefactors, although the dependence on parameters such as kf, I, and N c is 
correct. 



IV. EFFECTS OF TIMERE VERS AL- SYMMETRY BREAKING ON 

REFLECTION 

A. Change in < R > due to TR Symmetry 

Although individual diagrams for < Rij > cannot be evaluated precisely due to the 
inaccuracy of the diffuson approximation, we can evaluate the sum of all the diagrams over 
i,j from current conservation. Since we have proven the diagrams of Fig. 6 form a current- 
conserving set we may replace each diagram by its volume average after summing over 
In this case one is essentially using the relationship < R >= N c — < T > and calculating 
conductance diagrams. All diagrams except the one with the single crossed ladder vanish; 
this is the familiar diagram which upon evaluation in the cooperon approximation to the 
crossed ladder gives 5 < T >wl= — 1/3 |27] in the limit of L x 3> L y (quasi-ld sample). The 
only question raised by our derivation is whether one expects the cooperon approximation to 



be adequate here as again the points r, r' may coincide, and in fact standard treatments [28 
express the WL correction in terms of c(r, r). However the volume- averaging again appears 
to rescue this approximation. Any diagram with just a few scatterings will be of order 
l/L x , and may be neglected. It is only the diagrams with a number of impurity scatterings 
n ~ (L x /l) 2 which are long-range enough to give a result of order unity, and for such 
diagrams the vertex may indeed be treated as short range and the cooperon approximation 
is adequate. Since a weak magnetic field kills the cooperon contribution, the change in total 
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reflection due to breaking TR symmetry is 5 < R >=< R{B) > — < R(0) >= 
leads to a positive magnetoconductance which saturates at a value < g(B) > 
0) > =5g B = (+l/3)(e 2 /h). 



-1/3. This 
-<g(B = 



B. Evaluation of 5Ra 



Although we cannot evaluate all reflection diagrams precisely, we can accurately estimate 
the diagonal terms which correspond to the coherent backscattering effect discussed in the 
introduction. The simplest diagram is just the crossed ladder with no dressings and a fixed 
mode i,j at each vertex. For i ^ j this vanishes using Eq. (17) due to the orthogonality 
of different mode wavef unctions; for i = j the diagram is equal to the ladder sum for 
< Ru > except for the absence of the single-scattering diagram. Hence we have < 
R%% >coo P = 8Rd ~ 1- Thus we see that the diagonal reflection is enhanced in the presence 
of TR symmetry by roughly a factor of two, confirming our statement that this diagram is 
the mode-space manifestation of the coherent backscattering peak. This diagram makes no 
contribution to off-diagonal reflection. Therefore if coherent backscattering and WL were 
equivalent one would expect the WL effect in conductance to be e 2 /h and not (l/3)e 2 //i as 
we just found. The other dressed diagrams for < R t j > are not diagonal in mode-space and 
must add up to the difference (—2/3)e 2 /h. 

Unfortunately, due to the inaccuracy of the cooperon approximation there is no reason 
that explicit evaluation of these diagrams will yield precisely (— 2/3)e 2 /h. Thus we simply 
show that there exist diagrams with the appropriate order of magnitude and sign to give 
the expected effect. We note that all the other diagrams we have introduced which include 
dressings on the cooperon are of order 1/N 2 . Thus after summing over all modes they will 
indeed give a contribution of order unity. We evaluate the simplest of the diagrams (see 
Fig. 8) which yields a negative contribution to reflection coefficients. This will illustrate the 
general properties of these diagrams and also shows that time-reversal symmetry can indeed 
reduce the off-diagonal reflection (as well as enhancing the diagonal reflection). 

The diagram has an overall negative sign arising from the presence of two extra advance 
Green functions with no compensating retarded functions (each 1PGF gives a factor i). 
There are two such diagrams (the second is obtained by dressing the top of the diagram) 
with the same value. The cooperon is non-zero only if the transverse momentum on one 
side is the same magnitude as on the other side. This does not hold exactly for the crossed 
ladder, but it is still sharply peaked about reflection into the same mode. To lowest order 
in (kfl)^ 1 we have that the x-position of the impurity dressing must be less than that of 
both ends of the cooperon. The integration over the x-position of this point then will not 
introduce any additional phase. After integrating by parts and some algebra we find that 
this diagram and its conjugate contribute to reflection from mode i to mode j 




(22) 



After some simplifications this becomes 




(23) 
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As noted above, the contribution is negative and of order 1/N%. This diagram is typical 
of the off- diagonal contributions to reflection coefficients in the presence of time-reversal 
symmetry, i.e. it vanishes on volume averaging and is of order 1/N%, but gives a total 
contribution to reflection of the same order as the bare crossed ladder diagram because there 
is no restriction i — j. 

Such contributions are hard to understand from a (naive) semiclassical point of view. 
Since the diffuson and cooperon parts of such diagrams are typically represented in real- 
space, it has become common in the literature |]5| j20|j28[| to interpret each diagram as repre- 
senting a type of interference between paths. With this interpretation the cooperon repre- 
sents a time-reversed pair of counter-propagating paths which form a closed loop, consistent 
with the well-known interpretation of the CB peak mentioned in the introduction. However 
it must be kept in mind that the "paths" represented by such diagrams are not classical 
paths as they would be in a true semiclassical theory. Impurity-averaged perturbation the- 
ory in real- space allows any piece- wise linear path to contribute, and this complicates their 
semiclassical interpretation. For example, diagrams of the type shown in Fig. 6 in which 
the vertices are dressed with ladder insertions would correspond to two paths which diffuse 
together to the same point in the sample, but then diverge and traverse a loop in opposite 
directions, before possibly rejoining. Such a pair of paths is impossible classically, so the 
flux-sensitivity they impart to the reflection coefficients has no straightforward semiclassical 
interpretation. Although the white- noise model we have used has no simple classical limit, 
and one might speculate that it is in some respects untypical, it is noteworthy that such 
off-diagonal correlations are found numerically in simple chaotic systems ||||. Our work 
makes it clear that such contributions are essentially required by current conservation, and 
to our knowledge current conservation has not been proven for the semiclassical approxima- 
tion to the scattering coefficients. Thus we may hope that a deeper understanding of these 
results will arise from progress in the study of current conservation within the semiclassical 
approximation. 



V. SYMMETRY EFFECTS 



Our above analysis of the off-diagonal contribution to the WL effect in disordered con- 
ductors was motivated by numerical results obtained for the WL effect in ballistic conductors 
with chaotic classical dynamics These numerical studies also revealed another aspect 

of the WL effect which apparently has not been considered previously, the sensitivity of the 
effect to spatial symmetry. A particularly striking result of this work was a reduction of 
the WL positive magnetoconductance for ballistic junctions with double reflection symme- 
try; in this case the magnitude of the WL obtained numerically was slightly positive, but 
consistent with zero to the statistical accuracy of the calculations. Unlike the off-diagonal 
contribution to WL, which exists in both ordered and disordered media, the symmetry- 
induced contributions obviously do not exist in typical disordered media unless some very 
special process enforces a symmetry in the random potential. However it is still helpful 
to consider an ensemble of potentials which are random except for the presence of discrete 
spatial symmetries since in this case analytic calculations may be performed which provide 
insight into the origin of the sensitivity to symmetry and also may resolve the intriguing 
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question of whether certain symmetries can actually eliminate or reverse the sign of the WL 
effect. Thus we study the WL effect in symmetric "random" potentials in this section. 

It is possible to consider general spatial symmetry groups using the approach developed 
below. However, we find that only symmetry operations which are their own inverse give 
significant contributions; thus, we treat only single or double reflection symmetries. For 
definiteness we define the x-direction to be the direction of current flow and consider systems 
with reflection symmetry across the x-axis, the y-axis, or both. The symmetry operations 
on the random potential are represented by operators R x , R y and in the double reflection 
case R xy = R x R y . As before, we restrict ourselves to the two-dimensional case but extension 
to higher dimensions is trivial. The notation r^r i^. will be used to represent the point 
related to r by the relevant symmetry operator (the subscripts will be omitted when it is not 
necessary to distinguish the different symmetry operations). We modify the standard white- 
noise model for the impurity potential to include symmetry effects as follows: < V(r) >= 0, 
< V(r)V(r') >= CiU 2 [5(j — r') + X^f^r^ — r')], where the index /i runs over the complete 
set of symmetry operations. As before averages of higher even powers of V are given by the 
sum of all pairwise averages. 

With this modification the effective interaction mediated by the impurity potential no 
longer conserves momentum. If momentum p is lost by scattering off an impurity, p may be 
gained at the reflected impurity. As a result, the 1PGF is no longer diagonal in momentum 
space. However the modifications of the 1PGF will be negligible for our purposes. If we 
compare the contribution to the self-energy from scattering off the same impurity twice to 
that due to scattering from the impurity and its reflected impurity (see Fig. 9 ), the latter 
is only important if the impurity is within a few mean free paths of the reflection plane. 
This means that changes in the 1PGF produced by the new diagram will produce changes 
in the conductance that, after volume averaging, are lower order in l/L, where L is a typical 
sample dimension. Furthermore, even for two impurities lying near the symmetry plane one 
can verify that the phase of the Green function will not vary to first order in the transverse 
displacement of the impurity if the scattering is off the same impurity but will vary to first 
order when the scattering is from the reflected impurity (due to the length of the segment 
propagating from one to the other). As a result, the latter diagram will be reduced by 
the factor {k f l)- 1/2 . Thus we may use the usual impurity-averaged 1PGF in evaluating the 
conductance even in the presence of symmetry. 

Even though the 1PGF is unchanged the conductance is modified because new gener- 
alized ladder and crossed-ladder vertex corrections occur due to the presence of symmetry. 
An effect of this general type is fairly obvious from the semiclassical point of view in which 
the WL arises for constructive interference of paths with the same phase; in the presence of 
symmetry more such paths arise. In particular, if we interpret the real-space diagrams as 
representing paths, a crossed-ladder insertion in which the lower line is the precise reflection 
of the upper line (see Fig. 10) represents paths related by TR and reflection symmetry (or 
symmetries). We refer to this as a reflected crossed-ladder diagram and evaluate its con- 
tribution to the conductance in the cooperon approximation. These new diagrams arising 
due to symmetry in general have different numerical values than the standard WL diagram 
due both to changes in the value of the short-ranged current vertex and to changes in the 
value of the relevant cooperon insertion. We first evaluate the changes in the current ver- 
tex. We recall the standard procedure [EH] for evaluating the usual crossed-ladder diagram 
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(in the cooperon approximation) in real-space. The vertices of the diagram are treated as 
short-ranged and replaced with Vo5(r — r'), where Vq is the volume-average of the vertex. 
The diagram then becomes proportional to the integral of c(r, r) over the whole sample. To 
evaluate the coefficient V for the standard diagram, note that the vertices give a factor of 

G + (r, r 1 )^(r 1 )G"(r 1 , r')G"(r, r 2 H(r 2 )G + (r 2 , r') (24) 

After reversing the order of arguments to G (r, r 2 ) and G + (r 2 ,r') (which is permitted in 
the presence of time- reversal symmetry) and integrating over all the coordinates r, r', ri, r 2 
this becomes 

-Tr(G + v x G~G + v x G~). (25) 

In the reflected crossed-ladder diagram, one or more of the Green's functions in the 
vertex connect to a point r instead of r and similarly for r'. The new trace may then be 
written 

- Tr(G + v x G-R„G + v x G-R^, (26) 

where \i = x,y,xy. All the operators i? M commute with G since G(r, r') = G(r, r'); however 
while R x will also commute with v x , R y will anti- commute. In addition all the i? M satisfy 
R 2 — 1. Hence we may move one of the symmetry operators through the trace by com- 
mutation or anticommutation until it is adjacent to the other at which point the product 
becomes unity. Therefore we find that if /i = y,xy the coefficient of the new vertex is — V 
whereas if /i — x it is just V . 

The new conductance diagrams involve the reflected cooperon c(r, r). However it is 
easy to see that c = c because the kernel of the integral equation defining the reflected 
crossed-ladder G^ (r, r') ^ (r, is identical to that defining the usual crossed-ladder given 
the symmetry of the Green function under reflection. Hence the only difference from the 
standard diagram (other than the possible sign change of the vertex) is that the new diagrams 
are proportional to the volume average of c(r, r) rather than c(r, r). For a rectangular sample 
of length L x and width L y with the standard boundary conditions and no magnetic field 
we may write the cooperon in a spectral representation where the eigenf unctions of the 
diffusion equations are the normalized product of sine and cosine functions chosen to vanish 
at x = ±L x /2 and to have zero derivative at y = ±L y /2. Thus each term in the spectral 
sum will have the same value as in sum for c(r, r) or same magnitude but opposite sign, 
depending on whether the mode or modes reflected are even or odd under that reflection 
operation. For the most relevant case, that of a quasi-one-dimensional sample (L x ^> L y ), 
the cooperon is dominated by the terms with zero wavevector in the y-direction. In the 
standard cooperon sum these terms are proportional to 1/n 2 . If the symmetry present 

is either R y or R xy then this sum must be replaced by Y^=i(— l) n+1 /n 2 due to the alternating 
symmetry of the eigenfunctions in the x-direction. The first sum is equal to £(2) = 7r 2 /6, 
while the second is reduced by a factor of two. 

In addition to the fully-reflected crossed-ladder diagrams there are in principle other ver- 
tex corrections introduced by the imposition of symmetry. For example there are partially- 
reflected crossed-ladder diagrams in which some of the impurities on the lower Green function 
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line are reflected and some are not. It may be shown that these diagrams are lower order 
in (kfl) and l/L. Reflected ladder diagrams will give no contribution to the conductivity as 
can be seen by realizing that each vertex gives a contribution proportional to Tr(G + v x G~R) 
which clearly vanishes: if R anti-commutes with x, it trivially vanishes, while if R commutes 
with x, it vanishes after integrating over momenta. 

In summary, if R reflects across the y-axis or both the x-axis and the y-axis, for a quasi- 
1D sample, the diagram gives a contribution of opposite sign and one half the magnitude as 
the normal WL diagram. If R reflects across the x-axis only, the diagram contributes with 
the same magnitude and sign as the normal WL diagram. If there is 4-fold symmetry in the 
sample, the total WL effect is the sum of that for the non-symmetric case, plus diagrams in 
which R reflects across the x-axis alone, the y- axis alone, and both axes. As we reviewed 
above, the normal WL diagram is equal to —1/3 (in units of e 2 /h, neglecting spin); this 
gives the difference between the average conductance and the Drude average conductance 
(obtained from the diffuson diagrams). If we define Sg^, /i = x,y,xy as the difference 
between the Drude conductance and the average conductance in the presence of both TR 
and the relevant spatial symmetry our results imply: 

5g x = -2/3, 5g y = -1/6, Sg xy = -1/3. (27) 

Thus we arrive at several non-trivial conclusions. First, spatial symmetries do not reverse the 
sign of the WL effect; in all cases the conductance is reduced from its Drude value. Second, 
relative to the case of no symmetry, symmetry across the direction of current flow enhances 
WL, whereas symmetry perpendicular to the direction of current flow decreases WL. Finally, 
the combination of both symmetries leads to no change in the WL. All of these predictions 
can (and will) be tested numerically by averaging the conductance over realizations of the 
random potential which have the requisite symmetries and comparing them to the average 
without any symmetry. Thus we have in principle a means to "measure" the magnitude of 
the WL effect without imposition of a magnetic field. 

On the other hand in experiments on disordered conductors the WL effect is measured by 
the difference between the conductance at B=0 and at a magnetic field sufficient to destroy 
TR symmetry. If additional spatial symmetries are present the situation becomes more 
complex since in a uniform magnetic field a combination of spatial reflection plus TR may still 
yield an anti-unitary symmetry. This fact is familiar in the context of quantum chaos where 
it has been given the name "false TR symmetry breaking" ||29|| . In particular, a single spatial 
reflection (either R x or R y ) will reverse the sense of circulation(see Fig. 11) of any closed loop; 
since the TR operator will again reverse the sense of circulation (in a uniform field) a closed 
loop and its time-reversed, reflected partner will enclose the same flux and not be "dephased" 
by the field. Hence the single reflection diagrams will not be eliminated by a uniform 
magnetic field and will not contribute to the quantity 5gs =< g{B) > — < g(B = 0) >, 
which is the experimentally measured WL effect under normal conditions. The standard 
unreflected diagram is of course always eliminated by the field, so we conclude that the WL 
effect as usually measured will be unchanged by the imposition of single reflection symmetry 
and in quasi-lD will be Sgs = 1/3. However, double reflection symmetry will cause a 
closed loop to return to its original sense of circulation, hence TR plus double reflection 
will reverse the sense and such a diagram will be eliminated by a magnetic field. As shown 
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above, the double reflection diagram has one half the magnitude and opposite sign to the 
standard diagram, thus the measured WL effect 5gs = 1/6 and is reduced from its value 
in the absence of symmetry. Therefore the effect of the imposition of symmetry on the 
average conductance at zero field and on the field-dependent part of the conductance Sgs 
differ significantly. 

These results may be tested numerically in two contexts. First for disordered quasi- 
1D conductors with random potentials chosen to have no symmetry and then the sym- 
metries R x , R y , R xy . Since the average conductance with no symmetry already includes 
the standard WL effect it is more convenient to compare to the predictions for a quantity 
5g s =< g > sy mm ~ < 9 >non-symm- In addition we can compare the predictions of the 
theory to the numerical results obtained for ordered cavities with chaotic classical dynamics 
(the results which motivated this study). In this case we have no microscopic analytic theory 
at present. However, several studies pO"H32|l have used random matrix theory in either its 
Hamiltonian [BDJ or S-matrix [fni,B3] form to show that the standard (non-symmetric) WL 



effect is 0.25, and this is consistent with recent numerical results [p|,P|,pT|. Hence a plausible 
extension of our results for quasi-lD disordered systems would rescale 5g s , 5gs by the factor 
3/4; so for example we would expect a chaotic structure with double reflection symmetry to 
have Sgs = 1/8 instead of 1/6. In the chaotic case the averages are obtained by averaging 
over many values of the fermi energy as described in Ref. |J. Because of non-universal 
effects associated with short trajectories in the chaotic case, it is more difficult to compare 
two structures which are the same except for the imposition of spatial symmetry; however 
we have devised a method of analyzing the numerical results which we believe substan- 
tially overcomes this problem The numerical results were all obtained by the recursive 
Green function method |34] . |35|| which has been widely used in previous studies of disordered 
conductors and more recently in the study of ordered but chaotic microstructures ||. 

Our numerical results are presented in Tables I and II and compared to the analytic 
theory. It is worth noting several features of these results. First, the sign of the effect is 
always as expected from the theory; in particular the positive value of 5g s for symmetry 
perpendicular to the current as contrasted with the negative value for symmetry parallel 
is found both in the disordered and chaotic cases. Second, the magnitude of the effects 
found numerically tend to be smaller than that predicted by the theory, but are internally 
consistent. For example the theory predicts that Sgs should be the same for the three cases 
of no symmetry, symmetry parallel to current and symmetry perpendicular to the current. 
In both the disordered and chaotic systems these three cases closely agree even though the 
magnitudes found are 30 — 40% low |36| . Finally, the case of four-fold symmetry tends to 



have the greatest discrepancy from the theory, particularly in the chaotic structures. Overall, 
in our view the numerical results are in good qualitative agreement with the analytic theory 
and in reasonable quantitative agreement. 



VI. SUMMARY AND CONCLUSIONS 

We have presented a current-conserving approximation for the the local conductivity 
tensor which includes weak localization corrections and used it to show that the weak lo- 
calization correction to the conductance has important contributions not attributable to 
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coherent backscattering. The existence of such contributions, which do not have simple 
semiclassical analogues, complicates the interpretation of weak localization as compared to 
coherent backscattering. Numerical results indicate that such contributions occur both for 
disordered and chaotic conductors, raising questions about the possibility of a complete 
semiclassical theory of weak localization which are at present unanswered. A related ef- 
fect, relevant to weak localization in chaotic conductors, is the sensitivity of the effect to 
discrete symmetries. We showed that for disordered conductors reflection symmetries can 
either increase or decrease the weak localization effect, but cannot reverse its sign. Again 
these analytic results for disordered conductors are consistent with numerical results from 
chaotic conductors. The results relating to symmetry suggest the possibility of using the 
weak localization effect as a diagnostic of the symmetry of a ballistic microstructure. 
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APPENDIX A: EVALUATION OF REFLECTION COEFFICIENTS WITHIN THE 

DIFFUSON APPROXIMATION 

We evaluate the reflection coefficients within the diffuson approximation, using the 
boundary conditions of no net flux entering the sample, as discussed above. Since the 
purpose of this exercise is to demonstrate that the diffusion approximation is insufficient to 
handle reflection, we calculate only the total reflection coefficients to lowest order in (kfl)^ 1 , 
or, equivalently, to lowest order in 1/N C where N c is the number of transverse modes, equal 
to Lykf/n. The diagram which we must evaluate is shown in Fig. 7. The simplest method 
of evaluating this diagram is to use real-space in the x-direction and momentum-space in 
the y-direction as described in the text in section (3.4). We obtain in this case the same 
integral as in equation (|i~7D except that now we measure point x' from the same side of the 
sample as point x. In this case the diffuson, integrated over transverse positions, becomes: 



d(x, x') 



2 (x + 7il/4)(L x + nl/4-x') 

P Ly L x + nl/2 ' fOTX<X 

2 (x' + nl/4)(L x + nl/4-x) £ , (A1) 

p L " L x + vrZ/2 ' ^ x<x - 



To lowest order in (l/L x ) we may replace this by (2// 2 )(7r//4 + min(x, x')). For transmission 
coefficients we were guaranteed that one end of the diffuson would be to the right of the 
other due to the exponential decay of the Green's functions on the vertices; here we are not. 
The integration domain may be broken into x < x' and x > x' and the integral may be 
performed by parts. We obtain: 

^ 2 1 / 7T ViVj \ ... 

< Ra >=-— - + ^— . A2 

J iiN\4: (Vi + Vj)vf) V ; 
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Performing the sum over i,j numerically we obtain ps (0.73)7V c , instead of the value N c (up 
to corrections of order l/L) as required by current conservation and the known value of the 
average conductance. 

We know that calculating these diagrams using the correct ladder sum (instead of the 
diffuson approximation to it) will give the correct answer as we have proved that the ladder 
approximation conserves current. Since the exact ladder sum is not analytically tractable, 
we instead improve on the diffuson approximation by explicitly evaluating the single scatter- 
ing term in the ladder sum, and then evaluating a diagram with a single scattering dressing 
a diffuson. Thus the single-scattering contribution is evaluated exactly, and all double and 
higher scattering contributions are evaluated in the diffuson approximation. The single 
scattering can be done simply, giving a total contribution of ps 0.21iV c . The diffuson approx- 
imation for the second and higher scatterings is more tedious to evaluate but is found to sum 
to pa .63iV c . Hence the total average reflection in this improved approximation has increased 
to ^ 0.84iV c . In principle this method could be extended to even higher scatterings and an 
even better approximation for the average reflection will be obtained. 



20 



REFERENCES 



[1] E. Abrahams, P. W. Anderson, D. C. Licciardello, and T. V. Ramakrishnan, Phys. Rev. 
Lett. 42, 673 (1979). 

[2] For reviews of weak localization see G. Bergmann, Phys. Rep. 107, 11 (1984) and P. A. 

Lee and T. V. Ramakrishnan, Rev. Mod. Phys. 57, 287 (1985). 
[3] A. G. Aronov and Yu. V. Sharvin, Rev. Mod. Phys. 59, 755 (1987). 
[4] A. A. Abrikosov, L. P. Gorkov and I. E. Dzyaloshinski, Methods of Quantum Field 

Theory in Statistical Physics (Dover, New York, 1975). 
[5] S. Chakravarty and A. Schmid, Phys Rep 140, 195 (1986). 
[6] N. Argaman, unpublished. 

[7] A. Holle, J. Main, G. Wiebusch, H. Rottke, and K. H. Welge, Phys. Rev. Lett. 61, 161 
(1988). 

[8] H. U. Baranger, R. A. Jalabert, and A. D. Stone, Phys. Rev. Lett. 70, 3876 (1993). 

[9] H. U. Baranger, R. A. Jalabert, and A. D. Stone, Chaos, 3, 665 (1993). 
[10] K. Watson, J. Math. Phys., 10, 688 (1969). 
[11] D. A. de Wolf, IEEE Trans. Antennas Propag. 19, 254 (1971). 

[12] C.H. Lewenkopf and H.A. Weidenmuller, Ann. of Phys., 212, 53 (1991) and refs. therein. 

[13] M. P. van Albada and A. Lagendijk, Phys. Rev. Lett., 55, 2692 (1985). 

[14] P. E. Wolf and G. Maret, Phys. Rev. Lett., 55, 2696 (1985). 

[15] E. Akkermanns, P. E. Wolf and R. Maynard, Phys. Rev. Lett., 56, 1471 (1986). 

[16] B. W. Hapke, R. M. Nelson and W. D. Smythe, Science, 260, 509 (1993). 

[17] E. N. Economou and C. M. Soukoulis, Phys. Rev. Lett. 46, 618 (1981). 

[18] D. S. Fisher and P. A. Lee, Phys. Rev. B 23, 6851 (1981). 

[19] H. U. Baranger and A. D. Stone, Phys. Rev. B 40, 8169 (1989). 

[20] C. L. Kane, R. A. Serota, P. A. Lee, Phys. Rev. B 36, 5031 (1987). 

[21] C. L. Kane, P. A. Lee and D. P. DiVincenzo, Phys. Rev. B 38, 2995 (1988). 

[22] S. Hershfield, Annals of Physics 196, 12 (1989) and S. Hershfield and V. Ambegaokar, 

Phys. Rev. B 38, 7909 (1988). 
[23] M.V. van der Mark, M.P van Albada and A. Lagendijk, Phys. Rev. B 37, 3575 (1988). 
[24] E. Akkermans, P. E. Wolf, R. Maynard and G. Maret, J. Phys. France 49, 77 (1988). 
[25] A. Ishimaru, Wave Propagation and Scattering in Random Media (Academic Press, New 

York, 1978), Vol. 1, Chap. 9. 
[26] P. M. Morse and H. Feshbach, Methods of Theoretical Physics (McGraw-Hill, New York, 

1953), Part 2, Chap. 12. 
[27] P. A. Mello and A. D. Stone, Phys. Rev. B 44, 3559 (1991). 

[28] A. D. Stone, Physics of Nano structure, J. H. Davies and A. R. Long Eds., Institute of 

Physics Publishing, Bristol (1992) and references therein. 
[29] M. Robnik and M. V. Berry, J. Phys. A 19, 669 (1986). 

[30] S. Iida, H. A. Weidenmuller and J. A. Zuk, Phys. Rev. Lett. 64, 583 (1990) and Ann. 

of Phys. 200, 219 (1990). 
[31] H. U. Baranger and P. A. Mello, submitted to Phys. Rev. Lett., March 1994. 
[32] R. A. Jalabert, J.-L. Pichard, and C. W. J. Beenakker, submitted to Europhys. Lett., 

March 1994. 

[33] In the chaotic ballistic case, the numerical analysis was complicated by the fact that 



21 



changing the symmetry of the structure introduced a small but significant change in 
the classical transmission. In order to obtain results not influenced by the change in 
the classical transmission, we adopted the following procedure. (1) Fit a straight line to 
T(k). The slope, denoted T, has been previously shown to be the classical transmission 
probability ||. (2) Subtract the classical transmission, 5T(k) = T(k) — TkW/ir. (3) 
Average ST(k) over k for the different symmetry classes and compare them. 

[34] P. A. Lee and D. S. Fisher, Phys. Rev. Lett. 47, 882 (1981); D. J. Thouless and S. 
Kirkpatrick, J. Phys. C 14, 235 (1981); A. MacKinnon, Z. Phys. B 59, 385 (1985). 

[35] H. U. Baranger, D. P. DiVincenzo, R. A. Jalabert, and A. D. Stone, Phys. Rev. B 44, 
10637 (1991) and references therein. 

[36] One can think of many reasons that the agreement between numerics and analytic 
theory is imperfect by this amount for the chaotic conductors. For example, in this case 
our argument for neglecting the change in the 1PGF due to symmetry could fail as the 
size of the structure effectively sets the mean-free path. However it is puzzling that the 
results for the disordered case are consistently low by such an amount. There are certain 
parametric inequalities such as Z -C L x ,L y which are assumed in the analytic theory 
but are hard to satisfy strongly in the numerical simulations. Perhaps these account for 
the discrepancy; however in the case of the conductance fluctuations, which is similar 
in many respects, one finds agreement between analytic theory and numerics to 5% 
accuracy with little effort. 



22 



FIGURES 



FIG. 1. Self-consistent Born approximation (SCBA) to the average one particle Green's 
function. Dashed lines represent impurity scatterings, thin lines represent bare Green's functions, 
thick lines represent Green's functions after averaging. 

FIG. 2. Expansion for ladder sum. b) Expansion for crossed ladder sum. The + and - signs 
denote advanced and retarded Green's functions, L and X the ladder and crossed ladder sums. 

FIG. 3. Divergence of a general vertex in a conductivity diagram. The arrow represents the 
operation of taking the divergence, the wavy line at the vertex represents a current operator and 
the solid dot indicates a current vertex deleted by the divergence operation as described in the 
text. 

FIG. 4. a) Divergence of a vertex dressed with impurity scatterings, b) Divergence of bare 
vertex plus vertex dressed with ladder sum, represented by V. 

FIG. 5. The ladder approximation to conductivity. 

FIG. 6. Current conserving set of diagrams containing a crossed ladder with n scatterings. 
The triangular vertex denotes the dressed vertex V. a) The first set of diagrams as described in the 
text, b) The second set as described in the text. 

FIG. 7. Diffuson diagram used to calculate contributions to the transmission coefficients Tij 
to lowest order in 1/N C . 

FIG. 8. The simplest negative contribution to Rij due to time reversal symmetry, evaluated 
in the cooperon approximation. 

FIG. 9. The SCBA in the presence of symmetry. The corrections to the SCBA self-energy 
due to symmetry (last diagram) are found to be negligible. 

FIG. 10. Crossed ladder diagram with one side spatially reflected, denoted by the symbol X . 

FIG. 11. Effect of reflection on sense of a path. Single reflection reverses sense, double 
reflection preserves it. 

FIG. 12. Improved form of the diffuson approximation as discussed in the Appendix. 
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TABLES 



TABLE I. The effect of reflection symmetry on weak localization in quasi-lD disordered con- 
ductors. The results of numerical calculation are in good qualitative agreement with the analytic 



theory. The numerical results were obtained for the standard model [34,35| of a tight-binding 
Hamiltonian with on-site disorder. A rectangular strip with the following parameters was used: 
length/width = 8, width=120, disorder strength =0.9, kja = 1.61, mean-free-path=25, BA/4>o = 
or 10, and number of configurations=300. 

Symmetry Effects on Weak Localization: Quasi-lD Disordered Conductors 
Symmetry Type Sg s (theory) 5g s (num.) 5gs (theory) 5g% (num.) 

No Symmetry — 0.333 0.237 ± 0.026 

Symmetry |[ to Current (R x ) -0.333 -0.272 ± 0.037 0.333 0.211 ± 0.036 

Symmetry 1 to Current (R y ) +0.167 0.129 ± 0.037 0.333 0.273 ± 0.036 

4-fold Symmetry (R x R y ) 0.031 ± 0.051 0.167 0.059 ± 0.059 



TABLE II. The effect of reflection symmetry on weak localization in chaotic ballistic conduc- 
tors. The theoretical results are obtained by scaling the results for the disordered conductors by 
a factor of 3/4 while the numerical results are obtained for stadium- like billiards, the billiard in 
Fig. 1 of Ref. || for the asymmetric case and its symmetry generated analogs. For the numerical 
results, the average [33] was taken over the energy range kfW/ir € [4, 25] ([4, 15]) for the structures 



with (without) R x symmetry, the magnetic field used was BA/4>q = or 2, and kfa ph 1.0. Note 
the good qualitative agreement except in the four-fold symmetric case. 



Symmetry Effects on Weak Localization: Chaotic 


Ballistic Conductors 
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(num.) 


5g B (theory) 
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No Symmetry 









0.25 


0.147 
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0.018 


Symmetry to Current {Rx) 


-0.25 


-0.245 ± 


0.038 


0.25 


0.148 
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0.020 


Symmetry _L to Current {Ry 


+0.125 


0.170 ± 


0.030 


0.25 


0.161 
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0.019 


4-fold Symmetry (R x R y ) 





0.229 ± 


0.043 


0.125 


0.043 
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0.031 
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